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Abstract 



Measurements of the normalized rapidity (y) and transverse momentum (qj) dis- 
tributions of Drell-Yan muon and electron pairs in the Z-boson mass region (60 < 
Mu < 120 GeV) are reported. The results are obtained using a data sample of proton- 
proton collisions at a center-of-mass energy of 7 TeV, collected by the CMS experi- 
ment at the LHC, corresponding to an integrated luminosity of 36 pb _1 . The distribu- 
tions are measured over the ranges \y\ < 3.5 and qj < 600 GeV and compared with 
QCD calculations using recent parton distribution functions. Overall agreement is 
observed between the models and data for the rapidity distribution, while no single 
model describes the Z transverse-momentum distribution over the full range. 
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1 Introduction 

The production of Z and W bosons, which may be identified through their leptonic decays, 
is theoretically well described within the framework of the standard model. Total and differ- 
ential cross sections have been calculated to next-to-next-to-leading-order (NNLO) HUH]. The 
dominant uncertainties in the calculation arise from imperfect knowledge of the parton distri- 
bution functions (PDFs), from the uncertainty in the strong-interaction coupling oc s , and from 
the choice of quantum chromodynamics (QCD) renormalization and factorization scales. Mea- 
surements of the inclusive Z and W production cross sections performed by the Compact Muon 
Solenoid (CMS) experiment [3 J show agreement with the latest theoretical predictions both for 
the absolute value and for the ratios W + /W~ and W/Z. Likewise, agreement is found for the 
measurement of the dilepton mass distribution over a wide range [4j. 

In this paper, we present measurements of the normalized rapidity and transverse momentum 
distributions for Drell-Yan muon and electron pairs in the Z-boson mass region (60 < Mg£ < 
120 GeV). The results are obtained from a sample of proton-proton collisions at a center-of- 
mass energy of 7 TeV, recorded by the CMS detector at the Large Hadron Collider (LHC) in 
2010, and corresponding to an integrated luminosity of 35.9 ± 1.4 pb -1 . The measurement of 
the rapidity (y) and transverse momentum (qj) distributions of the Z boson provides new in- 
formation about the dynamics of proton collisions at high energies. The y distribution of Z 
bosons is sensitive to the PDFs, particularly when measured in the forward region (\y\ > 2.5), 
as done in this paper. The qj spectrum provides a better understanding of the underlying col- 
lision process at low transverse momentum, and tests NNLO perturbative QCD predictions at 
high transverse momentum. 

The rapidity is defined as y = \ In [(E + qi) I (E — qi)], where E is the energy of the Z-boson 
candidate and qi is its longitudinal momentum along the anti-clockwise beam axis (the z axis 
of the detector). The Z-boson y and qj are determined from the lepton momenta, which can be 
measured with high precision in the CMS detector. The measured differential dimuon and di- 
electron cross sections are normalized to the inclusive Z cross section, thereby canceling several 
sources of systematic uncertainties. 

The Z-boson y and qj distributions have been measured by the Tevatron experiments [jSHSI- 
In this paper, we report measurements which cover the range in rapidity up to 3.5 and in 
transverse momentum up to 600 GeV, a similar range to results recently reported by the ATLAS 
experiment |9j[l0|. The rapidity measurement is sensitive to the PDFs for proton momentum 
fractions (x) between 4 x 10~ 4 and 0.43. 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam- 
eter, providing a magnetic field of 3.8 T. Within the field volume are a silicon pixel and strip 
tracker, a crystal electromagnetic calorimeter (ECAL), and a brass/ scintillator hadron calori- 
meter (HCAL). The inner tracker measures charged particle trajectories in the pseudorapidity 
range \ t] | < 2.5 and provides a transverse momentum (pj) resolution of about 1-2% for charged 
particles with pj up to 100 GeV. The pseudorapidity rj is defined as t] = — ln(tan(0/2)), where 
9 is the polar angle with respect to the anti-clockwise beam direction. The electromagnetic cal- 
orimeter contains nearly 76000 lead-tungstate crystals that provide a coverage of \t]\ < 1.48 in 
a cylindrical barrel region and of 1.48 < < 3.0 in two endcap regions. The ECAL has an 
energy resolution of better than 0.5% for unconverted photons with transverse energies above 
100 GeV. The energy resolution is 3% or better for electrons with \tj\ < 2.5. The regions 
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3 Analysis procedure, data samples, and event selection 



(3.0 < \t]\ < 5.0) are covered by sampling Cherenkov calorimeters (HF) constructed with iron 
as the passive material and quartz fibers as the active material. The HF calorimeters have an 
energy resolution of about 10% for electron showers. Muons are detected in the range \rj \ < 2.4, 
with detection planes based on three technologies: drift tubes, cathode strip chambers, and re- 
sistive plate chambers. Matching segments from the muon system to tracks measured in the 
inner tracker results in a pj resolution of between 1 and 5% for muons with pj up to 1 TeV. Data 
are selected online using a two-level trigger system. The first level, consisting of custom hard- 
ware processors, selects events in less than 1 ^s, while the high-level trigger processor farm 
further decreases the event rate from around 100 kHz to about 300 Hz before data storage. A 
more detailed description of CMS can be found in Ref. |[TT|. 

3 Analysis procedure, data samples, and event selection 

The differential cross section is determined in each y or qj bin by subtracting from the number 
of detected events in a bin the estimated number of background events. The distributions 
are corrected for signal acceptance and efficiency and for the effects of detector resolution and 
electromagnetic final-state radiation (FSR) using an unfolding technique based on the inversion 
of a response matrix. The final result takes into account the bin width and is normalized by the 
measured total cross section. 

The measurements of the rapidity and transverse momentum spectra are based on samples of 
over 12000 Z-boson events reconstructed in each dilepton decay mode, and collected using 
high pj single-lepton triggers. The lepton identification requirements used in the analysis are 
the same as those employed in the measurement of the inclusive W and Z cross sections fl2|. 
For the Z-boson candidates selected, the pairs of leptons, I, are required to have a reconstructed 
invariant mass in the range 60 < M(? < 120 GeV. 

Muon events are collected using a trigger requiring a single muon, with a pj threshold that was 
increased from 9 to 15 GeV in response to increasing LHC luminosity during the data-taking 
period. The two muon candidates with the highest pj in the event are used to reconstruct a 
Z-boson candidate. Muons are required to have pj > 20 GeV, \r/\ < 2.1, and to satisfy the stan- 
dard CMS muon identification criteria described in Ref. [12]. In addition, the two muons are 
required to be isolated by calculating the sum of additional track momenta (I t rk) an d hadron cal- 
orimeter energy not associated with the muon (Ihcal) i n a cone AR = ^/{At]) 2 + (A(p) 2 < 0.3 
around the muon momentum direction, and requiring (J^ + Ihcal) /Pt(f) < 0.15. Infor- 
mation from the ECAL is not used as a criterion for isolation to avoid dependencies on FSR 
modeling [4J. At least one of the reconstructed muons must have triggered the event. The two 
muons in a pair are required to have opposite charges as determined by track curvature. The 
invariant mass distribution for selected events is shown in Fig. [I] We compare the kinematic 
distributions from the data and the simulations described below, and find that they agree. 

Electrons are detected in either the ECAL or the HF. For this analysis, the acceptance for elec- 
trons is defined to be within the fiducial region of ECAL, which overlaps with the silicon tracker 
region, or in the fiducial region of the HF. Electrons in this analysis can thus be observed over 
pseudorapidity ranges of \rj\ < 1.444 (ECAL barrel), 1.566 < \rj\ < 2.5 (ECAL endcaps), and 
3.1 < \t]\ < 4.6 (HF). The invariant mass distributions for selected events in the ECAL-ECAL 
and the ECAL-HF case are shown separately in Fig. [I] Events are selected online by a trigger 
requiring a single electron in the ECAL with pj > 17 GeV. The two electron candidates with 
highest pj in the event are used to reconstruct a Z-boson candidate, and at least one electron 
must be in the ECAL and have triggered the event. No requirement is applied on the charges 
of the electrons. Electrons are required to have pj > 20 GeV. Electrons reconstructed in the 
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Figure 1: Dilepton invariant mass distributions for the muon channel (top) and the electron 
channel (bottom). The bottom-left plot for the electron shows the invariant mass distribution 
for events with both electrons in the ECAL and the bottom-right for events with one electron in 
the ECAL and the other in the HE Each plot shows the data observation compared to the signal 
as predicted by POWHEG on top of the background estimated from a combination of simulation 
and data. The background is very low in the muon channel. 
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3 Analysis procedure, data samples, and event selection 



ECAL must have a matching track pointing to the reconstructed electromagnetic cluster and to 
be isolated and satisfy the general CMS electron identification criteria as described in Ref. fl2|. 

Electrons are reconstructed in the HF calorimeters from clusters of 3-by-3 towers centered on 
a seed tower with pj > 5 GeV. Each tower provides both a measurement of total energy de- 
posited and the energy deposited after the 12.5 radiation lengths (22 cm) of absorber closest 
to the interaction region. The two measurements are approximately equal for high-energy 
hadrons, while for electromagnetic particles the second measurement is typically a third of the 
total measured energy. Spurious signals from particles which pass directly through the pho- 
totube windows of the HF are rejected by requiring that the energy be shared among multiple 
towers. Electromagnetic clusters are selected by requiring the energy in the cluster to be at least 
94% of the energy in the 5-by-5-tower region containing the cluster. A further selection is per- 
formed using the ratio of the two energy measurements and the ratio of the two most energetic 
towers in the cluster to the total cluster energy. 

The detector acceptance is obtained from the simulation of the Drell-Yan process generated 
with the POWHEG US] matrix-element NLO generator interfaced with the PYTHIA (v. 6.422) HU 
parton-shower event generator, using the CT10 parametrization of the PDFs (TBI and the Z2 un- 
derlying event tune [TBI . The Z2 tune is the standard for CMS simulation and was tuned to the 
observed minimum-bias and underlying event characteristics at i/s = 7TeV [17]. The effect of 
FSR is simulated using pythia. In the muon channel, acceptance and efficiency calculations for 
the signal are performed using the full GEANT4-based |18| detector simulation, with additional 
smearing added to correct for observed differences in resolution between data and simulation 
as discussed below. For the electron channel, a parametrized simulation, matched to the reso- 
lution of the detector as measured in data, was used for efficiency and acceptance calculations. 
For the qj measurement, the electron acceptance is restricted to \rj\ < 2.1 to match the muon 
acceptance. 

The individual lepton detection and selection efficiencies are determined using a "tag-and- 
probe" method on the candidate lepton pairs. One of the leptons of the pair, the "tag", is 
required to pass all the selection requirements. The other lepton, the "probe", is selected with 
all requirements in the selection up to but excluding the requirement under study. The lepton 
pair is required to have an invariant mass consistent with the Z boson. When multiple tag- 
probe combinations are possible in a given event, one is chosen at random. The fraction of 
the probe leptons that also meet the requirement under study determines the efficiency of the 
requirement, after subtraction of the background from both samples using a fit to the dilepton 
invariant mass. In this manner, the efficiencies for the reconstruction, isolation, and trigger are 
measured sequentially. These efficiencies are compared with the efficiencies determined from 
the simulation to produce correction factors, some of which depend on the lepton kinematics, 
as discussed below. The efficiencies for an electron to form a cluster and a muon to form a basic 
track, both of which are very high, are taken from the GEANT4 simulation, which includes a 
modeling of inactive detector regions. The product of efficiency and acceptance for a given bin 
of y or qi is determined using Monte Carlo simulation as the ratio of the number of generated 
events reconstructed in the bin to the number of generated events before the FSR correction, 
using the single-lepton efficiencies determined from data. 

The single-muon trigger efficiency is determined separately for the different data-taking peri- 
ods and varies from 0.880 ± 0.008 at the beginning of the period to 0.924 ± 0.003 at the end. 
The single-muon trigger efficiencies are shown to be independent of px and tj within the accep- 
tance used in this analysis. The trigger efficiency for events with two muons of px > 20 GeV 
and \r]\ < 2.1 is 0.993 ± 0.005, averaged over data-taking periods. The average muon recon- 
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struction and identification efficiency for the selection used in this analysis is 0.950 ± 0.003. 
The uncertainty on the efficiencies is dominated by the data sample size for the tag-and-probe 
measurement. 

The single-electron trigger efficiency is measured to be between 0.96 ± 0.03 and 0.99 ± 0.01, 
varying as a function of pj and tj. For events with both electrons in the ECAL, the event 
trigger efficiency is greater than 0.999. For the electron channel, the total reconstruction and 
identification efficiencies determined from data range between 0.50 and 0.90 and are applied 
to the simulation as functions of pj and tj. Typical reconstruction and identification efficiency 
uncertainties are between 1 and 10%. However, the impact of these uncertainties on the final 
measurement uncertainty is greatly reduced by the normalization to the total cross section. 

The main sources of background in the measurement are Z — > tt and QCD multijet, tt, W + jets, 
and diboson production. Diboson production including a Z is considered to be a background 
for the measurement. All backgrounds except for QCD multijet production are evaluated us- 
ing Monte Carlo simulation. The Z — > tt events are generated with powheg+pythia. Events 
from tt, diboson production, and W + jets are generated using the MadGraph (v. 4.4.12) |(T9| 
matrix-element generator interfaced to pythia. Generated events are processed through the 
full GEANT4-based detector simulation, trigger emulation, and event reconstruction chain. We 
validated the use of the simulation to determine the background from the Z — » tt, tt, and dibo- 
son backgrounds by analyzing the qj spectrum for the ep pairs. These background processes 
are flavor-symmetric and produce twice as many ep pairs as ee or pp pairs. The analysis of this 
data sample matched the expectation from simulation. 

The QCD background is estimated using collision data samples. In the muon channel, the QCD 
background is estimated using a nonisolated dimuon sample corrected for the small contribu- 
tions in the nonisolated sample from prompt muons such as those from tt or Z-boson decay. 
The estimate is verified using a like-sign dimuon sample, since nonprompt sources of dimuons 
should have equal rates of like-sign and opposite-sign events. The QCD background in the 
muon channel is found to be very small. In the electron channel, the QCD background is larger 
and can be directly estimated by fitting the dielectron mass distributions in the data. The fit was 
performed in each measurement bin over the range 40 < M ee < 140 GeV using a linear combi- 
nation of a signal shape from simulation and a background shape determined by inverting the 
isolation and electron identification requirements in the data selection. 

After applying all analysis selection criteria, the total background fraction in the muon chan- 
nel is 0.4 ± 0.4%, consisting primarily of Z — > tt and tt processes and with an uncertainty 
dominated by statistical uncertainties in the background simulation. In the electron channel, 
the background fraction is 1.0 ± 0.5% for Z bosons reconstructed using two electrons in ECAL 
and 10 ± 4% for Z bosons reconstructed using one electron in ECAL and one in the HF, where 
the uncertainty is dominated by statistical uncertainties in the QCD estimate. In the electron 
channel, the QCD background is the dominant background component in every bin of rapidity 
and also at low qj. In the highest four qj bins, the Z — > tt and tt processes are the dominant 
contributions to the background. 

The bin width in rapidity (Ay = 0.1) is chosen to allow a comparison with previous measure- 
ments at lower center-of-mass energies. The bin widths in qj, which vary from 2.5 to 350 GeV, 
are chosen to provide sufficient resolution to observe the shape of the distribution, to limit bin 
migration, and to ensure a sufficient data sample in each measurement bin. 

The final measured y and qj distributions are corrected for bin-migration effects arising from 
the detector resolution and from FSR using a matrix-based unfolding procedure II20II . Large 
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4 Systematic uncertainties 



simulation samples are used to create the response matrices, which are inverted and used 
to unfold the measured distribution. This unfolding is applied to allow the combination of 
the muon and electron channels, which have different resolutions, and to allow the compar- 
ison with results from other experiments. The corrections resulting from detector resolution 
are calculated by comparing the generator-level dilepton distribution after FSR obtained from 
POWHEG+PYTHIA to that of the reconstructed simulated events, after smearing the momentum 
of each lepton with a parametrized function. The function is derived by comparing the Z mass 
distribution in data and the Monte Carlo detector simulation for different regions of rj and pj. 
In the muon channel, the smearing represents the observed difference between the resolution 
in data and in simulation. In the electron channel, as described above, a fully parametrized 
simulation is used for the acceptance and efficiency corrections. The corrections due to FSR are 
based on a Monte Carlo simulation using PYTHIA, and are obtained by comparing the dilepton 
y and qj distributions before and after FSR. These corrections are primarily important for the 
muon measurements, though large-angle FSR also has an impact on the electron distributions. 

4 Systematic uncertainties 

The leading sources of systematic uncertainty for the normalized distribution measurements 
are the background estimates, the trigger and identification efficiencies, the unfolding proce- 
dure, the calorimeter energy scale, and the tracker misalignment. Since the measurements are 
normalized by the total measured cross sections, several sources of systematic uncertainty can- 
cel, as they affect both the total rate and the differential rate in the same manner. For example, 
the uncertainty in the luminosity measurement cancels completely and uncertainties resulting 
from the lepton efficiencies and from the PDFs are significantly reduced. 

The muon and electron measurements share several theory-dependent uncertainties. Given 
the importance of the rapidity measurement in constraining PDFs, it is crucial to estimate the 
effect of the uncertainties in the PDFs on the determination of the bin-by-bin acceptance for the 
measurement. To evaluate the effect, we use the variations provided as part of the CT10 PDF 
set |Tl5| . For this PDF set, 52 variations are provided, each of which represents a shift in the 
PDFs by plus or minus one standard deviation along one of the 26 eigenvectors of the model. 
These eigenvectors are used to parametrize the uncertainties of the PDFs by diagonalizing 
the actual PDF model fit parameters, taking into account the unitary requirement and other 
constraints. The eigenvectors are not simply connected to specific observables, but represent 
an orthogonal basis in the PDF model space along which the uncertainties can be calculated. 
For each variation, the effect on the bin-by-bin acceptance normalized by the total acceptance 
was determined. The effects are combined in quadrature for each bin, separating negative 
and positive effects, to give the total uncertainty. The resulting uncertainties in the acceptance 
are less than 0.2% over the entire measurement range. However, the change in the shape of the 
distributions as a function of y is quite significant, up to 4% at high rapidity for some variations. 
These changes do not represent systematic uncertainties in the measurement - instead they 
represent the sensitivity of the analysis for constraining the PDFs. 

Several background processes, as described in Section |3j are predicted from Monte Carlo sim- 
ulation and compared with the data. A conservative estimate of the possible impact on the 
measurement is derived by varying the estimates of the small background from these sources 
by 100% based on the uncertainty due to the limited simulation sample size. We calculate the 
deviation of the central value of the normalized distribution in each bin when the background 
levels are varied. For the electron channel, the estimation of the bin-by-bin QCD background 
from data is a leading source of systematic uncertainty. Here, the error on the level of back- 
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Table 1: Fractional systematic uncertainty contributions for representative rapidity bins and 
transverse momentum bins in the electron and muon channels. 



\y\ Ran g e 


[0.0,0.1] 


[1.8,1.9] 


[3.0,3.1] 


Channel 


Muon 


Electron 


Muon 


Electron 


Electron 


Background Estimation 


0.002 


0.010 


0.002 


0.015 


0.047 


Efficiency Determination 


0.003 


0.005 


0.007 


0.007 


0.047 


Energy/ Momentum Scale 


0.001 


0.004 


0.001 


0.003 


0.009 


PDF Acceptance Determination 


0.001 


0.001 


0.001 


0.001 


0.001 


Total 


0.004 


0.012 


0.007 


0.017 


0.067 



Cjj Range 


[2.5GeV,5.0GeV] 


[110 GeV, 150 GeV] 


Channel 


Muon 


Electron 


Muon 


Electron 


Background Estimation 


0.004 


0.005 


0.019 


0.028 


Efficiency Determination 


0.010 


0.002 


0.010 


0.008 


Energy Scale 




0.022 




0.035 


Tracker Alignment 


0.015 


0.013 


0.023 


0.020 


Unfolding 


0.006 


0.004 


0.017 


0.001 


PDF Acceptance Determination 


0.002 


0.002 


0.001 


0.001 


Total 


0.020 


0.026 


0.036 


0.050 



ground in the signal region, 60 < M ee < 120 GeV, is dominated by the lack of data available in 
the dilepton invariant mass sideband regions. 

The trigger and the identification efficiencies are measured in the data as described above. The 
largest uncertainty in the efficiencies arises from the size of the data sample. To estimate the 
impact of these uncertainties on the final measurement, we change the efficiencies by plus or 
minus the amount of their statistical uncertainties and determine the changes of the normalized 
distribution. The changes from the central value are assigned as the systematic uncertainty aris- 
ing from the efficiency measurements, taking into account the cancellation effect from the rate 
normalization. The efficiencies from each stage of the selection are considered independently, 
and the resulting uncertainties are summed in quadrature. 

The systematic uncertainty from the unfolding procedure is estimated using alternative re- 
sponse matrices derived in several ways. We consider different generator models for the qj 
spectrum, which can affect the distribution of events within the bins. We vary the parameters of 
the detector resolution functions within their uncertainties. We also reweight the smeared spec- 
trum to match the data and evaluate the differences between the nominal and the reweighted 
unfolded spectra. In all cases, the effects amount to less than 0.5%. 

For the electron channel, the imperfect knowledge of the absolute and relative energy scales 
in the electromagnetic and forward calorimeters is a source of systematic uncertainty. Using 
Monte Carlo simulations, we estimate the effect of the scale uncertainties by scaling the ener- 
gies of electrons by amounts corresponding to the calibration uncertainties and the difference 
observed between the different calibration techniques used in the calorimeters. These energy 
scale uncertainties depend on the position of the electron within the calorimeters. We then 
determine the impact of these shifts on the observed distributions. 

The muon pj used in the analysis is based on the silicon tracker measurement. Thus, any mis- 
alignment of the tracker may directly affect the muon momentum resolution. The systematic 
uncertainty associated with tracker misalignment is calculated by reprocessing the Drell-Yan 
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5 Rapidity Distribution Results 



simulation using several models designed to reproduce the possible misalignments that may 
be present in the tracker. The bin-by-bin maximum deviation from the nominal Drell-Yan sim- 
ulation is used as estimator of the tracker misalignment uncertainty. In the electron channel, 
the sensitivity to the tracker alignment is determined by comparing the reconstructed y and pj 
using the calorimeter energy alone with those including the track measurements, for both data 
and simulation. 

The systematic uncertainties are summarized in Table [Tjfor representative values of y and qj in 
the muon and electron channels. After combining the effects discussed above, the total system- 
atic uncertainty in each bin is found to be significantly smaller than the statistical uncertainty. 

5 Rapidity Distribution Results 

The rapidity y of Z bosons produced in proton-proton collisions is related to the momentum 
fraction x + (x_) carried by the parton in the forward-going (backward-going) proton as de- 
scribed by the leading-order formula x± = £ ±!/ . Therefore, the rapidity distribution directly 
reflects the PDFs of the interacting partons. At the LHC, the rapidity distribution of Z bosons is 
expected to be symmetric around zero, therefore the appropriate measurement is the distribu- 
tion of Z bosons as a function of the absolute value of rapidity. The measurement is normalized 
to the total cross section (l/c da I d \y\), where a is the cross section determined by the sum of 
all observed y bins (|y| < 3.5), corrected to the total cross section as calculated from POWHEG 
with CT10 PDFs. The calculated correction between the measured and total y range is 0.983 
with an uncertainty of 0.001 from PDF variation. 

The measurements for the muon and electron channels are given in Table [2] and are in agree- 
ment with each other (reduced x 2 = 0.85) over the 20 bins where the measurements overlap. 
We combine these two measurements using the procedure defined in Ref. [21], which provides 
a full covariance matrix for the uncertainties. The uncertainties are considered to be uncorre- 
cted between the two analyses, since the only correlation between the channels is from the 
small PDF uncertainty. The combined measurements are shown in Table [2] and compared to 
the predictions made using CT10 PDFs in Fig.|2j 

To evaluate the sensitivity of this result to parameters of some of the more-recent PDF sets, 
we determine the change in the x 2 f° r each variation of the eigenvectors provided in the PDF 
sets. The CT10 PDF set has a x 2 of 18.5 for the base prediction, and the eigenvector-dependent 
changes in x 2 are shown in Fig. [3] The number of degrees of freedom (ndof) is 34. The 
MSTW2008 1)221 PDF set has a x 2 of 18.3 for its base prediction, and the eigenvector-dependent 
changes shown in Fig. |4] For both sets, several eigenvectors show significant sensitivity to our 
result, with CT10 showing a generally larger sensitivity. The HERAPDF 1.5 [23] PDF set, which 
has a x 2 of 18.4 for its base prediction, provides both eigenvectors and model dependencies as 
part of the PDF set. The changes in x 2 for both are shown in Fig. [5J The largest model depen- 
dencies with our measurement are the strange-quark PDF as a fraction of the down-quark-sea 
PDF. For the NNPDF 2.0 PDF set EH, the base prediction has a x 2 of 18.4. The NNPDF for- 
malism does not use eigenvectors, but rather replica PDFs sampled from the same space. In 
comparing our result with the 100 standard NNPDF 2.0 replicas, the majority have x 2 similar 
to the base, but some have x 2 values up to 34.5, indicating that these replicas are disfavored 
significantly by the new measurement. 



9 



Table 2: Measurement of the normalized differential cross section <jp J f° r Drell-Yan lepton 
pairs in the Z-boson mass region (60 < < 120 GeV) as a function of the absolute value 
of rapidity, separately for the muon and electron channels and combined. Detector geometry 
and trigger uniformity requirements limit the muon channel measurement to |y| < 2.0. The 
uncertainties shown are the combined statistical and systematic uncertainties. 





Normalized Differential Cross section 


V Ranee 


Muon 


Electron 


Combined 


[0.0,0.1] 


0.324±0.012 


0.359±0.015 


0.337±0.010 


[0.1,0.2] 


0.338±0.013 


0.326±0.016 


0.335±0.010 


[0.2,0.3] 


0.338±0.013 


0.344±0.017 


0.341±0.010 


[0.3,0.4] 


0.341±0.013 


0.355±0.017 


0.346±0.010 


[0.4,0.5] 


0.363±0.013 


0.339±0.017 


0.354±0.011 


[0.5,0.6] 


0.342±0.013 


0.351±0.018 


0.346±0.010 


[0.6,0.7] 


0.312±0.013 


0.360±0.018 


0.328±0.010 


[0.7,0.8] 


0.354±0.013 


0.331±0.018 


0.347±0.011 


[0.8,0.9] 


0.343±0.014 


0.355±0.018 


0.347±0.011 


[0.9,1.0] 


0.332±0.014 


0.332±0.018 


0.332±0.011 


1.0,1.1 


0.336±0.014 


0.316±0.018 


0.329±0.011 


[1.1,1.2] 


0.324±0.014 


0.352±0.019 


0.334±0.011 


[1.2,1.3] 


0.321±0.014 


0.332±0.019 


0.325±0.011 


[1.3,1.4] 


0.355±0.016 


0.321±0.019 


0.341±0.012 


1.4,1.5 


0.326±0.016 


0.313±0.019 


0.319±0.012 


[1.5,1.6] 


0.331±0.018 


0.330±0.020 


0.330±0.013 


[1.6,1.7] 


0.294±0.018 


0.306±0.022 


0.299±0.014 


[1.7,1.8] 


0.331±0.021 


0.332±0.024 


0.331±0.016 


[1.8,1.9] 


0.324±0.025 


0.294±0.024 


0.308±0.017 


[1.9,2.0] 


0.328±0.032 


0.328±0.026 


0.328±0.020 


[2.0,2.1] 




0.294±0.027 


0.294±0.027 


[2.1,2.2] 




0.298±0.029 


0.298±0.029 


[2.2,2.3] 




0.290±0.031 


0.290±0.031 


[2.3,2.4] 




0.278±0.035 


0.278±0.035 


[2.4,2.5] 




0.199±0.038 


0.199±0.038 


[2.5,2.6] 




0.249±0.040 


0.249±0.040 


[2.6,2.7] 




0.241±0.037 


0.241±0.037 


[2.7,2.8] 




0.256±0.035 


0.256±0.035 


[2.8,2.9] 




0.221±0.034 


0.221 ±0.034 


[2.9,3.0] 




0.165±0.035 


0.165±0.035 


[3.0,3.1] 




0.183±0.040 


0.183±0.040 


[3.1,3.2] 




0.228±0.045 


0.228±0.045 


[3.2,3.3] 




0.078±0.043 


0.078±0.043 


[3.3,3.4] 




0.105±0.051 


0.105±0.051 


[3.4,3.5] 




0.089±0.062 


0.089±0.062 
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5 Rapidity Distribution Results 
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Figure 2: The normalized differential cross section for Z bosons as a function of the absolute 
value of rapidity, combining the muon and electron channels. The error bars correspond to 
the experimental statistical and systematic uncertainties added in quadrature. The shaded area 
indicates the range of variation predicted by the POWHEG simulation for the uncertainties of 
the CT10 PDFs. 
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Figure 3: The change in % 2 when comparing the Z rapidity differential cross section measure- 
ment with the predictions of the NLO CT10 PDF set as each of the eigenvector input parameters 
is varied by plus or minus one standard deviation around its default value. 
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Figure 4: The change in % 2 when comparing the Z rapidity differential cross section measure- 
ment with the predictions of the NLO MST2008 PDF set as each of the eigenvector input pa- 
rameters is varied by ±90% confidence level (CL) around its default value. 
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Figure 5: The change in x 2 when comparing the Z rapidity differential cross section measure- 
ment with the predictions of the NLO HERAPDF 1.5 PDF set as each of the eigenvector input 
parameters (left) and the model parameters (right) is varied by one standard deviation around 
its default value. These together represent the full set of uncertainties in the HERAPDF 1.5 set. 
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6 Transverse Momentum Distribution Results 

Measurements of the qj distribution for Z bosons provide an important test of the QCD predic- 
tions of the initial-state gluon-radiation process. Perturbative QCD calculations are expected to 
provide a reliable prediction for the portion of the spectrum qj > 20 GeV, which is dominated 
by single hard-gluon emission. For qj < 10 GeV, the shape of the distribution is determined 
by multiple soft gluon radiation and nonperturbative effects. Such effects are simulated by 
Monte Carlo programs combining parton showering and parameterized models. These soft- 
gluon contributions can also be accounted for by resummation calculations in some Monte 
Carlo programs. 

For the qj measurement, the data are normalized to the cross section integrated over the accep- 
tance region \rj\ < 2.1 and pj > 20 GeV. The lepton pj and | rj | restrictions apply to both leptons 
of a dilepton pair. The restriction on the electron pseudorapidity (compared to that used for the 
rapidity measurement) allows the combination of the two channels and a more straightforward 
physics interpretation, as the two measurements refer to the same rapidity range and have the 
same PDF dependence. 

Table 3: Measurement of the normalized differential cross section for Drell-Yan lepton pairs in 
the Z-boson mass region (60 < M(j < 120 GeV) as a function of qj, separately for muon and 
electron channels and for the combination of the two channels. The distribution is normalized 
by the cross section for Z bosons with both leptons having \t]\ < 2.1 and pj > 20 GeV. The 
uncertainties listed in the table are the combined statistical and systematic uncertainties. 



qj Range ( GeV) 


Muon Channel 


Electron Channel 


Combination 


[0.0,2.5] 


(3.21 ±0.14) 


x 10" 


-2 


(3.24 ± 0.25) 


x 


io- 


-2 


(3.22 ± 0.13) 


X 


10" 


-2 


[2.5,5.0] 


(5.89 ±0.21) 


x 10" 


-2 


(6.03 ± 0.32) 


X 


io- 


-2 


(5.92 ± 0.17) 


X 


io- 


-2 


[5.0,7.5] 


(5.51 ±0.20) 


x 10" 


-2 


(5.32 ± 0.32) 


X 


io- 


-2 


(5.50 ± 0.16) 


X 


io- 


-2 


[7.5,10.0] 


(3.90 ±0.18) 


x 10- 


-2 


(4.20 ± 0.30) 


X 


io- 


-2 


(3.96 ± 0.14) 


X 


io- 


-2 


[10.0,12.5] 


(3.49 ±0.16) 


x 10- 


-2 


(3.60 ± 0.28) 


X 


io- 


-2 


(3.53 ± 0.12) 


X 


io- 


-2 


[12.5,15.0] 


(2.74 ±0.15) 


x 10- 


-2 


(2.70 ± 0.25) 


X 


io- 


-2 


(2.72 ± 0.12) 


X 


io- 


-2 


[15.0,17.5] 


(2.23 ±0.14) 


x 10- 


-2 


(2.00 ± 0.22) 


X 


io- 


-2 


(2.16 ±0.10) 


X 


io- 


-2 


[17.5,20.0] 


(1.68 ±0.12) 


x 10- 


-2 


(1.59 ±0.20) 


X 


io- 


-2 


(1.65 ± 0.09) 


X 


io- 


-2 


[20.0,30.0] 


(1.14 ±0.04) 


x 10- 


-2 


(1.20 ±0.05) 


X 


io- 


-2 


(1.16 ±0.04) 


X 


io- 


-2 


[30.0,40.0] 


(6.32 ± 0.28) 


x 10- 


-3 


(5.62 ±0.31) 


X 


io- 


-3 


(5.98 ± 0.27) 


X 


io- 


-3 


[40.0,50.0] 


(3.53 ±0.21) 


x 10- 


-3 


(3.18 ±0.24) 


X 


io- 


-3 


(3.38 ± 0.18) 


X 


io- 


-3 


[50.0,70.0] 


(1.74 ±0.10) 


x 10- 


-3 


(1.90 ±0.12) 


X 


io- 


-3 


(1.81 ±0.09) 


X 


io- 


-3 


[70.0,90.0] 


(7.76 ±0.71) 


x 10- 


-4 


(7.86 ± 0.77) 


X 


io- 


-4 


(7.79 ± 0.54) 


X 


io- 


-4 


[90.0,110.0] 


(4.87 ±0.55) 


x 10- 


-4 


(4.57 ± 0.59) 


X 


io- 


-4 


(4.75 ± 0.42) 


X 


io- 


-4 


[110.0,150.0] 


(1.79 ±0.22) 


x 10- 


-4 


(2.18 ±0.26) 


X 


io- 


-4 


(1.93 ±0.17) 


X 


io- 


-4 


[150.0, 190.0] 


(7.10 ± 1.40) 


x 10- 


-5 


(4.82 ±1.31) 


X 


io- 


-5 


(6.00 ± 0.99) 


X 


io- 


-5 


[190.0,250.0] 


(1.17 ±0.51) 


x 10- 


-5 


(2.05 ± 0.64) 


X 


io- 


-5 


(1.51 ± 0.43) 


X 


io- 


-5 


[250.0,600.0] 


(2.24 ± 0.78) 


x 10 


-6 


(0.81 ±0.52) 


X 


io- 


-6 


(1.29 ± 0.44) 


X 


io- 


-6 



The measurements from the muon and electron channels are tabulated in Tableland are found 
to be compatible with each other over the full qj range (reduced x 2 = 0.74). The combina- 
tion of the muon and electron results is also performed following Ref. [21J. The alignment 
uncertainty is treated as correlated between the two channels, and other uncertainties are 
treated as uncorrelated. The combined measurement is presented in Fig. |6j where the data 
points are positioned at the center-of-gravity of the bins, based on the POWHEG prediction. For 
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Figure 6: The Z-boson transverse momentum distribution found from combining the muon 
and electron channels, compared to the predictions of the POWHEG generator interfaced with 
pythi A using the Z2 tune. The error bars correspond to the statistical and systematic uncertain- 
ties added in quadrature. The band around the theoretical prediction includes the uncertainties 
due to scale variations and PDFs. The horizontal error bars indicate the bin boundaries, and the 
data points are positioned at the center-of-gravity of the bins, based on the POWHEG prediction. 
The inset figure shows the low qj region on a linear scale. 
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qi > 20 GeV, we compare the data and the prediction of POWHEG+PYTHIA with the Z2 tune 
and find ^ 2 /ndof = 19.1/9, where ndof is equal to the number of points minus one because 
of the normalization. We have taken the full covariance matrix into account when computing 
the x 1 values. At low momentum, there is poor agreement, suggesting the need for additional 
tuning of the combination of POWHEG and PYTHIA in this region, where both contribute to the 
observed qj. 

At low transverse momenta, i.e. qj < 30 GeV, the distribution is determined by nonperturba- 
tive QCD, which is modeled by PYTHIA with a few free parameters. Several parameter sets 
called "tunes" are available, including the Perugia 2011 |25], ProQ20 11261 , and Z2 tune IIT6I . 
The shapes predicted with these tunes are compared to this measurement in Fig. [7] Agree- 
ment is observed for the Z2 (# 2 /ndof = 9.4/8) and the ProQ20 tunes (;\; 2 /ndof = 13.3/8), 
but disagreement for the Perugia 2011 tune (^ 2 /ndof = 48.8/8) and for POWHEG+PYTHIA 
(^; 2 /ndof = 76.3/8). These results provide a validation of the Z2 tune for a high momentum- 
scale process that is rather different from the low-momentum-scale processes that determine 
the characteristics of minimum-bias events and the underlying event from which the parame- 
ters of the Z2 tune were originally obtained. 

At high qj, the precision of the prediction is dominated by the perturbative order of the calcu- 
lation and the handling of the factorization and renormalization scale dependence. In Fig.|8]the 
measured normalized differential distribution is compared to the prediction of POWHEG as well 
as the "Fully Exclusive W, ZProduction" (fewz) package J27) for q T > 20 GeV and \t]\ < 2.1, 
calculated at both 0(a s ) and 0(a 2 ). The predictions were each normalized to their own pre- 
dicted total cross sections. The FEWZ calculation used the effective dynamic scale definition 

J M| + (^t) 2 rather than the fixed scale of the Z-boson mass. The fewz 0(a 2 ) prediction pro- 
duces a ^ 2 /ndof of 30.5/9, which is a poorer agreement than the POWHEG prediction (19.1/9), 
particularly at the highest qj. 

7 Summary 

Measurements of the normalized differential cross sections for Drell-Yan muon and electron 
pairs in the Z-boson mass region (60 < Mfj < 120 GeV) have been reported as functions of 
dilepton rapidity and transverse momentum separately. The results were obtained using a 
data sample collected by the CMS experiment at the LHC at a center-of-mass energy of 7 TeV 
corresponding to an integrated luminosity of 36 pb -1 . The rapidity measurement is compared 
with the predictions of several of the most recent PDF models and the agreement evaluated as 
a function of the PDF set eigenvectors. An overall agreement between the models and the data 
is observed. The measured transverse momentum distribution is compared to various tunes 
of the PYTHIA generator for low transverse momentum and to 0(a s ) and 0(a 2 ) predictions for 
high qj. No single model describes the normalized differential cross section of the Z transverse 
momentum over the full range. These measurements significantly extend previous Tevatron 
results and complement recent LHC results in rapidity and transverse momentum. 

Acknowledgments 

We wish to congratulate our colleagues in the CERN accelerator departments for the excellent 
performance of the LHC machine. We thank the technical and administrative staff at CERN and 
other CMS institutes. This work was supported by the Austrian Federal Ministry of Science and 
Research; the Belgium Fonds de la Recherche Scientifique, and Fonds voor Wetenschappelijk 



16 



7 Summary 



^ 0.08 

i 

% 0.07 
0^0.06 
-D 0.05 
^ 0.04 



LI 



0.03 P^H 
0.02 

0.01F- 





~i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



L dt = 36 pb" 1 at \/s = 7 TeV 



CMS : 



L_. 



• data (e + ji combined) 

POWHEG+PYTHIA Z2 

PYTHIA Z2 

PYTHIA ProQ20 



pythia Perugia 2011 



|<2.1, p >20 GeV 








20 25 

Qj [GeV] 



30 



Figure 7: The combined electron and muon measurement of the Z-boson transverse momentum 
distribution (points) and the predictions of several PYTHIA tunes and POWHEG interfaced with 
PYTHIA using the Z2 tune (histograms). The error bars on the points represent the sum of 
the statistical and systematic uncertainties on the data. The lower portion of the figure shows 
the difference between the data and the simulation predictions divided by the uncertainty S 
on the data. The green (inner) and yellow (outer) bands are the ±15 and ±2$ experimental 
uncertainties. 



17 




20 30 40 50 100 200 300 



Q T [GeV] 

Figure 8: The combined electron and muon Z-boson normalized differential cross section 
as a function of transverse momentum (points) and the POWHEG and FEWZ predictions for 
qj > 20 GeV. The horizontal error bars indicate the bin boundaries and the data points are po- 
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upper plot represent the uncertainty on the predictions from factorization and renormalization 
scales and PDFs. The lower plot shows the ratio between the data and the theory predictions. 
The bands in the lower plot represent the one standard deviation combined theoretical and 
experimental uncertainties. 
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